Nonlinear coupling of nano mechanical resonators to Josephson quantum circuits 
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We propose a technique to couple the position operator of a nano mechanical resonator to a 
SQUID device by modulating its magnetic flux bias. By tuning the magnetic field properly, either 
linear or quadratic couplings can be realized, with a discretely adjustable coupling strength. This 
provides a way to realize coherent nonlinear effects in a nano mechanical resonator by coupling it 
to a Josephson quantum circuit. As an example, we show how squeezing of the nano mechanical 
resonator state can be realized with this technique. We also propose a simple method to measure the 
uncertainty in the position of the nano mechanical resonator without quantum state tomography. 



Introduction. Historically, mechanical systems have 
not been the favorite proving ground of quantum mechan- 
ics because of their macroscopic nature. This situation 
has changed in recent years thanks to the impressive ad- 
vance of nano fabrication technologies. It is now possible 
to make nano mechanical resonators with frequencies of 
gigahertz Q] and quality factors approaching 10 5 at milli 
Kelvin temperatures [2|. This opens the possibility of 
studying coherent quantum behavior in mesoscopic me- 
chanical systems One can also exploit the quantum 
properties of mechanical degrees of freedom for applica- 
tions in areas such as weak force detection |||, precision 
measurement 0, and quantum information processing 
0- Recently, evidence for quantized displacement in a 
nano mechanical resonator has been observed || . 

All micro and nano mechanical devices require some 
means of transduction. As discussed in 3], an excel- 
lent way to engineer and detect the quantum modes of a 
nano mechanical resonator is to couple the resonator to 
Josephson device based solid state circuitry on which 
quantum coherent control has been demonstrated. A 
straightforward coupling scheme for this purpose is to 
voltage bias the resonator and use the position dependent 
electrostatic interaction between the nano resonator and 
a charge island. This is discussed in [j| and used in all 
previous studies (see 01 f° r an incomplete list) . Though 
conceptually simple and practically feasible, this scheme 
has the limitation that the dominant term in the induced 
coupling is always linear in the resonator position. Con- 
sequently, it can only be used to realize linear effects in 
the nano resonator. As is well known in quantum op- 
tics 0], nonlinear effects are indispensable for thorough 
study and control of the dynamics of harmonic oscilla- 
tors. They are required to produce essential processes 
such as squeezing and parametric amplification. They 
are also known to be necessary for universal quantum in- 
formation processing on continuous variables |l2| . Thus, 
it is highly desirable to develop new coupling schemes 
which can introduce nonlinear effects into the nano res- 
onator system. 

In this work, we propose an attractive alternative 
method to couple the nano mechanical resonator to a 
Josephson quantum circuit by modulating the flux bias of 



a SQUID device. Our scheme makes it possible to realize 
both linear and nonlinear processes on the nano mechan- 
ical resonator. We show how we can generate squeezing 
of the nano mechanical resonator using this method and 
propose a simple way to measure the reduction in the 
uncertainty in the resonator's position. 

The device and its working principle. Our scheme is 
illustrated in Fig. ^ ( a )- Here, the nano mechanical 
resonator is in one arm of a SQUID device. The SQUID 
loop is biased with a perpendicular magnetic field B, such 
that the flux threading the SQUID is the product of B 
and the area of the SQUID loop. 
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FIG. 1: (a) A nano mechanical resonator in one arm of a 
SQUID device. The SQUID loop is threaded by a magnetic 
field. Its flux bias will then depend on the position of the 
nano mechanical resonator, (b) A more complicated design 
in which a nano mechanical resonator is coupled to a charge 
qubit. The shaded box is the superconducting island of the 
charge qubit. 

We see that the area of the SQUID loop is dependent 
on the position of the nano mechanical resonator. In Fig. 
H(a), we denote the width and length of the SQUID loop 
W and L. The length and small displacement of the nano 
mechanical resonator are I and X, defined such that the 
area of the SQUID loop is WL + IX. The total flux bias 
of the SQUID loop is $° + BIX, where $° = BWL is 
the flux bias corresponding to the equilibrium position 
of the nano mechanical resonator. With the phase drops 
of the two junctions being <j>i and (f>2, the Josephson en- 



ergy of the SQUID is Ej = -E 7 cos 



-2E°j cos( 



i - E°j cos e> 2 = 
i - 2 )/2cos((/>i + (j) 2 ), where E°j = I c ^ /2tt 
is the Josephson energy of the (identical) junctions, I c 
the junction critical current and $o the flux quantum. 
Since the superconducting order parameter is singly val- 
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ued, we must have fa - fa + 27r($" + BlX)/<& = 2-np 
for some integer p leading to 



Ej = -2E°jC0s{ 



7T$° ttBI 



X) cos< 



(1) 



where = (fa + fa)/ 2 is the average phase across the 
junctions. 

It is clear from Eq. |Q that in general the Joscphson 
energy of the SQUID is a nonlinear function of X, the 
position of the nano mechanical resonator. Treating X 
as a small perturbation, we find that there are two flux 
bias points of particular interest: = (n + l/2)$o an d 
<I>p = n$o, where n is an integer. When = (n + 
l/2)$0) to lowest order in X the Josephson energy of 
the SQUID is Ej = (-l) n 2Ej(ir(n + l/2)l/WL)Xcos(l>. 
In this case, Ej has a linear dependence on X and thus 
linear coupling between the nano resonator and SQUID 
is realized. On the other hand, when the SQUID loop is 
biased at $® = n$ , to lowest order in X the Josephson 
energy is 



Ej 



-(-l) n 2Ej{l - (-nnl/WL) 2 X 2 /2} cos0. (2) 



Here, the Josephson energy of the SQUID has a quadratic 
dependence on the position of the nano mechanical res- 
onator. 

From the perspective of micro and nano mechanical 
engineering, our coupling scheme can be considered a 
magnetic transduction method. Compared to the elec- 
trostatic transduction scheme the distinctive advan- 
tage of our scheme is that the coupling can be either 
linear or quadratic, depending on the flux bias point of 
the SQUID. Also, notice that the coupling strength can 
be adjusted (discretely), since Ej depends on n as in Eq. 
(J2J. It is then possible in principle to operate in both 
linear and nonlinear, as well as both weak and strong, 
coupling regimes. In practice, the SQUID will be part 
of a Josephson quantum circuit (for instance a charge 
or flux qubit), and the modulation scheme above then 
couples the nano mechanical oscillator to the Josephson 
quantum circuit. 

Note the nano mechanical resonator must be supercon- 
ducting for our scheme to work. This can be realized by 
using a metalized resonator 0, 0] as long as the mag- 
netic field does not exceed the critical field strength of 
the superconductor. 

Squeezing of the nano mechanical resonator. An 
important nonlinear effect on a harmonic oscillator is 
squeezing. It is the suppression of the uncertainty in 
the position (or momentum) of the oscillator at the price 
of increased uncertainty in the conjugate variable. For 
nano mechanical resonators, reduction in the position 
uncertainty is important in applications such as weak 
force detection yj. Also, some coherent nonlinear pro- 
cess such as squeezing is essential for universal quantum 
information processing on continuous variables 12]. Our 



nonlinear coupling scheme makes such coherent nonlin- 
ear processes inaccessible by previous methods 0, 0] 
possible. This provides a way of introducing nonlinear- 
effect-induced operators, for instance, as gates within a 
complicated quantum circuit. Previously, generation of 
certain squeezed states of nano mechanical resonators has 
been studied by a few authors 0>0- However, since 
no coherent nonlinear transduction scheme on the nano 
mechanical resonator was known, these schemes had to 
function incoherently. They used dissipation and mea- 
surement to generate the needed nonlinearity and were 
incapable of effecting unitary evolutions. 

In the following we show how we can realize squeez- 
ing on a nano mechanical resonator using our scheme to 
couple it to a charge qubit ^t|- First, we notice that the 
simple prototype circuit in Fig. ^( a )j though convenient 
in illustrating the essence of our scheme, has a few dis- 
advantages. In particular, according to Eq. J2J), when a 
nonlinear coupling is realized the Josephson self-energy 
is also the largest (2Ej), therefore we do not have com- 
pletely independent control over the coupling strength 
and Josephson self-energy. To overcome this and other 
disadvantages, we consider the design shown in Fig. 
(b). Here, we have two identical SQUIDs (feft and right) 
biased at equal but opposite fluxes <& e / = — <I> er = BWL 
(this can be realized by twisting the arms of one of the 
SQUIDs, for instance). The nano mechanical resonator is 
in the arm of one of the SQUIDs and coupled to its phase. 
The big loop is biased with an external flux & x . This 
structure bears some similarity to the scalable charge 
qubit scheme in 0]> an d the advantage is that both the 
flux biases of the individual SQUIDs and the big loop 
can be tuned. 

In order to realize nonlinear coupling between the 
nano mechanical resonator and the charge qubit, we 
bias the SQUIDs at $ e ; = -$ er = ri<E> . Us- 
ing the same argument that leads to Eq. (J2J), we 
can derive the Josephson energy of the circuit in 
Fig- G] (b). To lowest order in the resonator 
position X, Ej = -(-1)"4E^ cos^^o) cos + 
(— l) n Ej (n nl / W L) 2 X 2 cos 4> i, where 4>i and fa are the 
phases of the left and right SQUIDs (the averages of 
the phases of the two junctions in the SQUIDs), and 
<t> = (4>i + fa)/ 2 is the average phase of the SQUIDs 
conjugate to the charge number on the island. If 
we bias the big loop at <& x — (2m + l/2)$o, m an 
integer, only the coupling term survives, so Ej = 
(~l) n E j (irnl /WLfX 2 cos cj>i. 

The charge island possesses an adjustable gate volt- 
age V g which is applied through the gate capacitance 
C g . When it is biased close to n g — C g V g /2e = 1/2, 
the states with and 1 excess Cooper pairs comprise 
the low energy Hilbert space of the qubit. Considerin 
these charge states the spin up and spin down states 
in an effective two state system, we can use the Pauli 
matrices to describe operators acting on the system, and 
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its uncoupled Hamiltonian takes the form E z o~ z /2, where 
E z = (2n g — l)(2e) 2 /2Ct and Ct is the total capacitance 
of the charge island. 

As usual, the nano resonator is treated as a harmonic 
oscillator with position operator X — 5Xo(a + af) Q, 
where a is the annihilation operator, 8X0 = y/h/2MuJo 
is the zero point fluctuation in the resonator's position 
X, and M and luo are the mass and frequency of the 
resonator. Its uncoupled Hamiltonian is ujoa'a. 

The system Hamiltonian is then H = E z o- z /2+u!oa^a— 
(X n /2)(a + a)) 2 (jy, where the coupling strength A„ = 
-(~l) n E j(TTnlSX /WL) 2 . If we choose E z = 2uj and 
shift to the rotating frame defined by E z a z /2 + Lu^a 'a, 
we obtain the following Hamiltonian in the rotating frame 

El: 

H R = i{X n /2){a 2 a+ - a^cr') + H„ , (3) 

where = (<j x ±ia y ) /2 and are off resonance terms 
of magnitude X n /2 oscillating at frequencies of 2loq and 
higher. Since the realizable coupling strength X n is usu- 
ally much smaller than the resonator frequency ujq, we 
can adopt the rotating wave approximation to drop i? Wo 

El. 

In addition to the operating point discussed above, we 
also consider another set of bias conditions in which the 
SQUIDs are biased at flux and E z is changed (by tun- 
ing the gate voltage V g ) to E' z = 2ujq + 5E Z . We also 
bias the big loop slightly away from (2m + 1 /2)<&o us- 
ing a small ac field, Q x = (2m + 1 /2)<&o + S§ x cos2woi 
where it8§ x /&q <C 1. In this case the charge qubit 
is decoupled from the resonator. To lowest order in 
ttS$ x /<5>o, the system Hamiltonian in the same rotat- 
ing frame is H — 5E z a z /2 + E x a x /2 + H' Uo , where 
E x = 4E°j(ir5& x /&o) and the rapidly oscillating term 
H' Uo will have negligible effect if E x is small compared 
to 4wn and the system evolves for an appropriate dura- 
tion |l9j. In the above Hamiltonian both 8E Z and E x of 
the qubit can be adjusted 0]; therefore we can perform 
arbitrary rotations on the state of the charge qubit. In 
particular, if we choose large values for 8E Z and E x , we 
can realize pulsed operations on the charge qubit and flip 
its state quickly. 

We now consider a spin echo like process in which we 
let the system evolve under the Hamiltonian J2J for short 
periods of time of duration At. In between each such 
time interval we apply a quick ir pulse a x to the charge 
qubit to flip its state, so that the evolution for two pe- 
riods is governed by exp{— iH^At}a x exp{— iH^Ai}a x . 
Since a x exp{(A„At/4)[(a 2 — a^ 2 )a x + i(a 2 + a) 2 )a y ]}a x 
= exp{(A„At/4)[(a 2 - a^ 2 )a x - i(a 2 + a^)a y }}, 
this evolution operator can be simplified: 
exp{— iHjiAt}a x exp{— iHjjAt}a x « exp{(A„At/2)(a 2 — 
a^ 2 )o~ x }. If we initialize the charge qubit in the a x — 1 
state, and repeat this procedure N times, the evolution 



operator on the state of the resonator becomes 

5( K ) = exp{|(a 2 - a t2 )}5 (4) 

where n = X n NAt. This is a squeezing operator on 
the nano mechanical resonator with squeezing param- 
eter k. Under the squeezing operator, a transforms 
to S* (k)ciS(k) = a cosh k — sinh(K) and it can be 
shown that the position uncertainty decreases exponen- 
tially El : Ax = V( x2 ) - ( x ) 2 = AA(0)e~ K . 

In the above, we used a multi-loop circuit which allows 
bias and control of the individual loops. Such structure 
and control are easy to realize and widely used in cur- 
rent Josephson quantum circuit design and experiment 
El EH- Also, high-precision spin echo control of 
superconducting qubits has been realized experimentally 
|22j . Sophisticated microwave pulse sequences can be ap- 
plied and it was observed that the decoherence time of 
the superconducting qubit increases due to the spin echo 
control. Therefore, our scheme is within the reach of 
current technologies. 

Ideally, one would like a large coupling strength in or- 
der to generate appreciable squeezing in a short period 
of time. For a magnetic field B « 0.2T, I « 30^m, 
5Xq sa 5 x 10~ 13 m, and a critical current of about 60nA, 
the coupling strength A„ is about 4MHz. 

Decoherence. Unlike previous methods [TH , our 
scheme effects a coherent quantum process. The influ- 
ence of decoherence must be considered in detail. For 
this purpose, we use the Master equation El 

dp/dt = -i[H sq , p]+£(a, 7„, N n )+£(cr~ , 7,, N q )+C(a z , 74,, N q ), 

where p is the density matrix of the nano resonator - 
charge qubit system, H sq — iX(a 2 — a' r2 )a x /2 is the ef- 
fective squeezing Hamiltonian, j n — loq/Q is the decay 
rate of the nano mechanical resonator determined by its 
quality factor Q, 7 g and 7^ are the relaxation and de- 
phasing rate of the charge qubit, N n and N q are the 
mean values of the bath quanta dependent on temper- 
ature 11], and the Liouvillian operator C(A,j,N) = 
{l/2)-f(N+l){2ApA^-A^ Ap-pA^ A)+(l/2)-fN(2A^ pA- 
AA^ p — pAA^). Among the various decoherence sources, 
the dephasing of the charge qubit is dominant. For most 
nano mechanical resonator frequencies 0, we only bias 
the charge qubit slightly away from the charge degen- 
eracy point (\2n g — 1| <C 0.1) and the dephasing time 
T 2 is greater than 100ns Using the Master equa- 

tion we can derive equations for the expectation values 
of the dynamic variables of the system which are then 
solved. In Fig. |3 (a), we plot the time dependence of 
the nano resonator position uncertainty AX. We have 
used the following conservative set of experimental pa- 
rameters: resonator frequency loq/2it = 250MHz, quality 
factor Q = 10 4 , temperature T = 20mK, squeezing pa- 
rameter A = 5MHz, with 7^ and 7^ chosen such that the 
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relaxation time T\ — 1/is and dephasing time T2 — 100ns 
|23j| . Initially the nano resonator is assumed to be in a 
thermal equilibrium state and the charge qubit is in the 
a x = 1 state. 

It is clear from Fig. (a) that appreciable squeezing 
can be generated even when the charge qubit 's dephasing 
rate is severe (twice the squeezing parameter) , indicating 
that the squeezing is robust against decoherence. As time 
progresses, the squeezing becomes less effective and even- 
tually AX starts to increase. This is easy to understand 
from the effective Hamiltonian i\(a 2 — a <f2 )(7 x leading to 
Q; as the charge qubit dephases, (a x ) decreases and the 
squeezing effect weakens and eventually disappears. Not 
surprisingly, the maximum squeezing achievable increases 
with decreasing dephasing rate, as shown in Fig. 0(b). 

(a) (b) 
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FIG. 2: (a) Time dependence of the nano resonator position 
uncertainty AX (normalized to the initial uncertainty) under 
the effect of squeezing and decoherence. See text for experi- 
mental parameters and initial conditions, (b) The maximum 
achievable squeezing (minimum AX) as a function of the de- 
phasing rate. All other parameters are the same as in (a). 

Measurement. Once we generate squeezing on the nano 
mechanical resonator, a scheme is needed to measure the 
uncertainty in its position to confirm the squeezing effect. 
One way to do this is quantum tomography which mea- 
sures the resonator's Wigner function. As shown in [25) . 
the Wigner function can be measured by displacing the 
state of the harmonic oscillator, letting it interact with a 
two state system, and measuring the polarization of the 
two state system. Though this method can determine 
all the information about the oscillator, its direct appli- 
cation in our solid state system is hindered by technical 
difficulties. To displace the resonator's state arbitrarily, 
we will need couplings to both the position and momen- 
tum operators of the resonator with continuously variable 
relative strengths, which is not easy to realize in the nano 
resonator. Also, in order to calculate AX, we need the 
Wigner function over its entire parameter space making 
it necessary to sweep through a large parameter space. 
Large displacement will take a long time to effect and the 
decoherence will affect the measurement result. 

Here we propose a simplified method to measure the 
mean and variance of the position operator of the nano 
resonator. It is based on the measurement of the gener- 
ating function Tr{pe tKX ). We notice that (X) and (X 2 ) 



can be calculated from the generating function by 

(X) = Tr(pX) = -iA-Tr(pe iKX )\ K=0 (6) 

and 

(X 2 ) = Tr(pX 2 ) = -^Tr(pe^ x )\ K=0 . (7) 

Therefore, measuring the generating function in the 
vicinity of k — allows us to determine the position 
uncertainty of the nano mechanical resonator AX 2 . 

The generating function can be measured by first 
preparing the charge qubit in the o~ z = 1 state. Then 
turning on a strong coupling XXa x between the charge 
qubit and the resonator, for a time t, will cause the 
charge qubit states corresponding to a x = ±1 to acquire 
a phase shift ^fXtX. Next, the interaction is turned off 
and the gate H exp{i(7r/4)(n-er)} is applied to the charge 
qubit, where H = (a x + a z )/^/2 is the Hadamard gate, 
n = (0, — cos 77, sin 77), r\ an angle chosen to be or n/2. 
We then measure the polarization of the charge qubit, 
(a_z) = P<r z =i — P a :,=-!, which can be shown to equal 

E3 

Tr{pRe{e"'e lKX }), (8) 

where k = 2Xt. Choosing rj = and ir/2 then yields the 
real and imaginary part of the generating function, which 
in turn allows us to calculate AX 2 by Eqs. © and J7J. 

This method based on the measurement of the generat- 
ing function has a few attractive advantages. It requires 
only strong linear coupling of the charge qubit to the 
position operator of the resonator, which as discussed 
before can be easily realized in our scheme by biasing 
the SQUIDs at (n + l/2)$ and the big loop at flux. 
We only need to measure in the vicinity of k = 0, there- 
fore the measurement can be done quickly, implying less 
influence by decoherence. Since the polarization of the 
charge qubit can be measured with high fidelity [24| , our 
scheme is realistic given currently available technology. 

Conclusion. We have proposed a scheme to couple a 
nano mechanical resonator to Josephson quantum cir- 
cuits by modulating the magnetic bias of a SQUID. This 
allows us to realize coherent nonlinear effects on the nano 
resonator, which are essential for but so far missing in 
the study of nano mechanical resonators. Though we fo- 
cused on the squeezing of a nano mechanical resonator 
by coupling it to a charge qubit, our scheme can be eas- 
ily tailored to other purposes and adapted for coupling 
to other Josephson quantum circuits. It can be directly 
extended for quantum manipulation of multiple nano res- 
onators. Then, entanglement can be generated and con- 
firmed using the simple measurement scheme. This can 
provide a practically feasible approach for unambiguous 
demonstration of quantum behavior in nano mechanical 
resonators. 
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